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Abstract. We consider alternative theories of gravity with a direct coupling 
between matter and the Ricci scalar We study the relation between these theories 
and ordinary scalar-tensor gravity, or scalar-tensor theories which include non- 
standard couplings between the scalar and matter. We then analyze the motion 
of matter in such theories, its implications for the Equivalence Principle, and the 
recent claim that they can alleviate the dark matter problem in galaxies. 



1. Introduction 

An explanation for the nature of the unknown form of energy, usually dubbed dark 
energy, which appears to constitute 76% of the energy budget of the universe 1 and is 
responsible for its late time accelerated expansion is still pending. The dominant role 
of gravity in the dynamics of the late time universe and the unexpected characteristics 
of dark energy, namely the fact that it violates the Strong Energy Condition and may 
even be as exotic as phantom energy, have lead some to the idea that dark energy 
might be nothing but an indication of the inability of General Relativity (GR) to 
describe gravity consistently at cosmological scales According to this approach, 
some other theory of gravity, which actually constitutes an extension of GR, might 
correctly describe the cosmological dynamics without dark energy and, at the same 
time, reduce to the usual GR phenomenology at scales where the latter has already 
been very successful. 

Several theories have been put forth, often as toy models with the purpose 
of providing proofs of principle that such an explanation to the recent puzzles of 
cosmology can indeed be feasible. Most of them also seem to draw some motivation 
from high energy physics, as they are presented as candidates for a low energy effective 
action of some more fundamental theory {e.g., String Theory). A typical example is 
that of the so-called f{R) gravity [3]. These theories, which are proposed as minimal 
extensions of GR, are described by the action 



where g is the determinant of the metric gab, n = SirG, f{R) is a general function 
of the Ricci scalar R, L„i is the matter action, and ip collectively denotes the matter 
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fields (we adopt the notations and conventions of Ref. [1]). When f{R) is taken to 
be linear, this action reduces to the usual Einstein-Hilbert action. Then, Einstein's 
equations can be derived with standard metric variation, but also with an independent 
variation of the metric and of the connection [31 U] called Palatini variation, under the 
crucial assumption that the matter action does not depend on the, now independent, 
connection [5]. Nevertheless, when f{R) is not linear in i?, the two variational 
principles lead to different theories: metric and Palatini f{R) gravity [J respectively 
(see O [6l [71 m [9] for early works and reviews) . Additionally, one could decide to allow 
the independent connection to enter the matter action when the Palatini variation is 
used (which appears to be in accordance with the geometrical meaning of a connection 
to define the covariant derivative [5]), in which case the outcome is yet another distinct 
version of f{R) gravity: metric-affine f{R) gravity [l0| . 

Much attention has been paid to f{R) theories of gravity lately, especially 
concerning the cosmological evolution but also the various viability constraints [TT] 
(with metric-affine f{R) gravity being probably the exception). An extension of metric 
f{R) gravity has been recently proposed in Ref. [12], in which the Ricci scalar R 
acquires an explicit coupling to the matter Lagrangian. The action of the theory is 

j d\^g\^^ + [l + \h{R)]L„}^ , (2) 

where is the matter Lagrangian, /i 2 are (a priori arbitrary) functions of the Ricci 
curvature i?, and the usual coefficient of R in the Einstein-Hilbert action (SttG, where 
G is Newton's constant) has been absorbed in fi{R). 

This action, and variants of it, were considered for various purposes: in Ref. |12| . 
([2]) was studied as an alternative to dark matter because it is claimed to give rise to 
phenomenology similar to MOND gravity [l^ on galactic scales; in [Ml [7], a variation 
of ([2]) was studied as an alternative to dark energy by setting /i (i?) = R and keeping 
only the nonminimal coupling of matter to the Ricci curvature, in order to explain the 
current acceleration of the universe. Along the same line of thought, the idea of making 
the kinetic term of a (minimally coupled) scalar field dependent on the curvature, 
while keeping fi (R) = R was exploited in attempts to cure the cosmological constant 
problem [TSllIS] (see also [T71 [T5] V In light of its relation with a scalar-tensor theory 
explained below, the action ^ is reminiscent also of the one describing a nonminimally 
coupled electromagnetic field [T^ . 

Variation of the action ([2]) with respect to the metric yields 

Fi{R)Rab - \h{R)9ab - V,VbFi(i?) + gabOF^iR) 



= -2XF2{R)L„,Rab + 2A(VaVfc - gabD)L,nF2{R) 



+ [1 + Xf2{R)]Tab , (3) 

where Fi 2{R) = /{ 2(^)7 the prime denotes differentiation with respect to the 
argument, and 

rr . = _^ S{V-gL,^) 

^ab ■ ^ > 



X In the Palatini formalism, the connection T^j^^ is independent of the metric. Therefore, R is replaced 
in the action |(TJ by 7^ = g°'^'R-abi where TZab is the Ricci tensor constructed with this independent 
connection. 
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As a consequence of the explicit coupling of the matter Lagrangian density to the 
Ricci curvature, there is an energy exchange between matter and gravity beyond the 
usual one always present in curved spaces, which is reflected in the non-vanishing 
of the covariant divergence of the matter stress-energy tensor Tab- The corrected 
conservation equation assumes the form [121 El E] 

^''Tab = , , ' {L,ngab - Tab) V' R . (5) 
t + A/2 

The fact that the stress energy tensor is manifestly not divergence-free (in the 
representation used to write down the action ([2])) can be interpreted as a violation 
of the so-called metric postulates fIU\ . This seems to imply that the theory exhibits 
violations of the Equivalence Principle (EP). This, as well as the fact that one could 
potentially tune the parameter A to reduce the effects of such violation below current 
experimental accuracy, has already been mentioned in Ref. jT2]. Nevertheless, to the 
best of our knowledge, no detailed study of this feature of the theories described by ([2|) 
has been carried out so far. Even though the metric postulates can be a good criterion 
for constructing theories which do not violate the EP, one can not safely conclude that 
any theory that does not satisfy the metric postulates in some particular representation 
indeed violates the EP. First of all, the metric postulates are representation-dependent 
statements, i.e., not invariant under field redefinitions j21l . In general, there might be 
other representations of the theory in which the metric postulates are indeed satisfied 
More importantly, even in theories that indeed exhibit violations of the EP, such 
as those considered here, the metric postulates or the non-conservation of energy do 
not themselves provide quantitative estimates of the deviations from the EP. 

The purpose of this paper is exactly to shed some light on these issues. We 
begin by examining, in section [21 the relation between the theories described by 
the action ordinary scalar-tensor gravity, and scalar-tensor theories with an 
anomalous coupling between matter and the scalar field. It is shown that the former 
can not be cast into the form of any of the latter (as usually done in "ordinary" 
f{R) gravity without the i?-matter coupling). This not only imphes that the theory 
under consideration can not be cast away as a known theory, but also that the R- 
matter coupling can not be eliminated with simple field redefinitions. We proceed, in 
section |3l to examine the effect of the i?-matter coupling on the motion of particles 
and fields. It is shown that massless particle trajectories and field propagation are 
actually unaffected (at least at high frequencies). More remarkably, even the motion 
of massive objects (such as perfect fluids or test particles) seems to remain unaffected. 
This implies that detecting deviations from geodesic motion in the theories under 
investigation might be more difficult than expected. On the other hand, it also casts 
doubts on whether they can actually account for the phenomenology for which they 
were introduced in |12j . 

2. Is i?-coupled modified gravity a deja vu? 

It is customary in metric (and also Palatini) /(i?) gravity to introduce auxiliary scalar 
fields, in order to re- write the action as a scalar-tensor theory [24] . Let us check what 

§ A characteristic example is scalar-tensor gravity: the metric postulates are satisfied in the Jordan 
frame but not in the Einstein frame. However, there are no violations of the Einstein EP, since the 
latter is a characteristic of the theory, not of its representations 1211 (see 1221 1231 for a discussion 
of the Einstein EP in nonminimally coupled scalar field theory, which is a particular realization of 
scalar-tensor gravity). It is just that in one frame this is made manifest. 
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would happen if we followed the same procedure here: we begin by introducing two 
scalars ^' and 4> and considering the action 

S'^ j d^x^g\^^ + [\ + \h{4>)]L„,, + ^{R-4>)^ . (6) 

The variation of this action with respect to 4" yields R = (f> and one recovers the 
action Now, varying ([6]) with respect to <j) yields 

* = + A/^i™ . (7) 

Assuming that at least one of /1.2 is non-linear in R (/" ^ or f2 7^ 0, or both) [jj] 
and replacing eq. ^ back in eq. ([6]) one gets (see also [l4l[7]) 

S'= J d^x ^ + [1 + A/2(0)] 

i/{ + A/^L„,) (i?-0)| . (8) 

This action is obviously dynamically equivalent to the action In the special case 
A = 0, in which ([2]) reduces to the action ([T]), and ([8]) reduces to 

S' = j d^x^g\^^+L„, + \f[{R-<P)Y (9) 

or, with the simple field redefinition $ = /((0) and the introduction of the potential 
y($) = {(j)f'i — /i)/2, it can be rewritten as a Brans-Dickc theory with a potential 
and vanishing Brans-Dicke parameter wq |24j : 

S'^ j d^x ^ - V{^) +L„Y (10) 

However, when A 7^ 0, one can see directly from the action ([8]) that the coupling 
between matter and R docs not disappear. Even if wc perform the field redefinition 
mentioned earlier and introduce the functions V"(<i>), J7($) = 1 + A(/2 — 0/2), and 
X($) = /2, at best we can write 

S' ^ J d^xV^[^~V{<f) + U{<S>)L,n + XX{<i>)RL„,} ■ (11) 

This is clearly not an ordinary scalar-tensor theory: first, there exists an unusual 
coupling between the scalar and the matter in the third term. This is reminiscent 
of extensions of scalar-tensor gravity which include similar couplings, such as the 
theories considered by Damour and Polyakov [26^. However, the presence of the R- 
matter (-scalar) coupling in the last term distinguishes the action pT|) even from these 
generalized scalar-tensor theories. Alternative ways to introduce a single scalar are 
also unable to alleviate the i?-matter coupling (see the Appendix). 

Also, one should not be mislead to think that the i?-matter coupling can be 
eliminated, judging from action ((6]). Even though in eq. ([6|) there is indeed no explicit 
coupling between R and the matter, there is a coupling between 5* and R. '5 in 
turn is just an auxiliary scalar which is algebraically related to the matter fields 
through eq. Consequently, the i?-matter coupling is just made implicit through 

II The particular case in which both functions are linear in R, f" = /j' = 0, is burdened with serious 
viability issues |25| and will not be studied further. 
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the introduction of an auxiliary (non-dynamical) field. Note the difference with cf) (or 
$) which does carry dynamics (it is algebraically related to R but not to the matter). 
Therefore, we conclude the following: 

• Introducing scalar fields helps in avoiding the presence of non-linear functions of 
R in the action (as in ordinary f{R) gravity). 

• However, the theories under scrutiny cannot be written into the form of a scalar- 
tensor theory with a minimal coupling to matterl^ I.e. the i?-matter coupling 
cannot be eliminated, or replaced by some unusual couplings between matter and 
the scalar field (it is, therefore, not an artifact of some peculiar representation). 

• One could "hide" the i?-matter coupling by considering the theory as a multi- 
scalar tensor theory (starting from action ^ and possibly considering also field 
redefinitions for the two scalars). However, in this case special attention should 
be paid to the fact that one of these scalar does not actually carry dynamics. 
This, for instance, distinguishes this theory from the multi-scalar-tensor theories 
already considered in the literature [27j . 

Consequently, the theories described by the action ^ can not be cast into the form 
of a scalar- or multi-scalar-tensor gravity which has already been extensively studied 
in the literature and its phenomenology is well known: The presence of the i?-matter 
coupling which is exactly the subject of interest here, is what distinguishes these 
theories and deserves further attention. 

3. i?-coupling, geodesies, and EP violations 

As already mentioned, the explicit coupling to the Ricci curvature described by the 
action ([2]) can potentially lead to non-geodesic motion and violations of the EP, judging 
from the fact that the stress-energy tensor is not divergence-free, as expressed by 
eq. (O. Although this is pointed out in Ref. [T^] with the caveat that such violations 
could be controlled by the value of the parameter A, this issue has not yet been 
addressed in detail. In this section we perform a more thorough analysis. First, we 
consider massless particles and high frequency fields, and then we move on to perfect 
fluids. 

3.1. Massless fields 

In [12], it is stated that the explicit coupling of the matter Lagrangian to the Ricci 
curvature does not alter the null geodesic equation: however, this statement is not 
proved, and the reader might be left with the impression that it is motivated by the 
resemblance of the action ^ with the one of Einstein frame scalar-tensor theories, 
which we have shown to be fallacious. As a matter of fact, the authors of [12] are more 
interested in the correction to the worldlines of massive particles and the corresponding 
MOND-like phenomenology; however, this aspect of the theory should not be left 
unchecked. 

^ By scalar-tensor theory we mean a theory with a single scalar field mediating gravity. For theories 
with more than one scalars we use the term multi-scalr-tensor theory. 
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3.1.1. Null dust: The equation for null geodesies can be derived in a straightforward 
manner from the conservation equation of a null dust fluid (see e.g.^ Ref. [H]). 
Therefore, if the i?-coupling were to induce any corrections to the null geodesic 
equation, these would show up in this derivation. 

The stress-energy tensor of a null dust is 

Tab^PUaUb, UcU^^O, (12) 

where p is the fluid energy density and u'^ is the four- velocity of masslcss fluid particles. 
Since the perfect fluid is an "averaged" and not an exact description for matter, it is 
more common in this case to work directly with the stress-energy tensor and avoid 
any reference to a Lagrangian. Unfortunately, in our case this is not possible since 
Lm enters explicitly in the conservation equations ((S]). Even though more than one 
choices for the Lagrangian have been used in the literature [301 El] > the most natural 
(and general) choice seems to be simply the pressure Lm = P [30] The corrected 
conservation equation ^ then becomes 

\F 

pu'^VbUa + pUaV''Ub + UaUbV'' p = - pUau'^V bR . (13) 

t + A72 

Now this equation can be written in the form 

U^\JbUa = Qua , (14) 

where 

1 XF 

e = -v^fc + -MbVV-^ — ^u''VbR- (15) 
p 1 + A/2 

Eq. ([T^ states that the four-velocity is parallel-transported along the path: this is the 
deflnition of a geodesic curve. The fact that this equation differs from the more familiar 
form u'^Vcu" = is simply because the latter is affinely parametrized. Therefore, the 
explicit i?-coupling does not change the equation of null geodesies. 

3.1.2. Masslcss scalar field: Let us consider now a massless scalar field (j) described 
by the Lagrangian density and energy-momentum tensor |32j 

L,„ = -iv,0V^0, (16) 

Tab = Va-/) VbC^ - ^ 5afcVe</> V^(/. . (17) 

With a correction proportional to gabLm — Tab = ^^a4'^b't'j (non-)conservation 
equation ([5]) is 

(V^V,0) Vf,0 + (Va0) (VaV^0) (V,0) - i (V^c/.) Va^c^ 

^-—^[yaC^){ybC^)y'R. (18) 
t + A/2 

By projecting onto V°0, one obtains the corrected Klein-Gordon equation 

+ T— (^'"^) V'i? = . (19) 

t + A/2 

+ Note that the choice of the Lagrangian is particularly meaningful here: in GR these different choices 
would anyway lead to the same field equations, unlike here where the actual expression for Lm enters 
the field equation of the theory explicitly. Therefore, all results strongly depend on this choice. 
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By taking the high frequency hmit 

(/.(x") = (/.o(x'^)e'^(-^) , (20) 

with the phase S a rapidly varying function of x"^ and the amphtude 0o i^"") slowly 
varying, and neglecting the gradients and second derivatives of (/>o, eq. (|19p implies 

- SaS" + iV'Sc = - tt4t ^Sb^'R , (21) 

^ + A/2 

where Sc = VcS* is the phase gradient and the tangent to the worldline of the scalar 
particle in the geometric optics approximation. This equation yields 

SaS" = (22) 

{i.e., the spacetime trajectory is null), and 

V^^, = - Tt4t ^b^'R ' (23) 

which expresses the fact that the "scalar photon" is not transversal unless A = or 
is orthogonal to the gradient of R. By further applying the covariant derivative 
operator to eq. and using the fact that [Vq, V^] 5 = 0, one obtains 

S'^WaS'' = . (24) 

Therefore, the worldlines of scalar particles in the high frequency approximation are 
null geodesies. This conclusion would not change if one were to consider instead a 
massless, non-minimally coupled, scalar field described by the Lagrangian density 

im = -^V,0V^0-|i?02 , (25) 

where ^ is a dimensionless coupling constant. In fact, the corresponding Klein-Gordon 
equation 

XF 

t + Aj2 

contains the extra term —£,R4> that can be interpreted as a tidal effect on the scalar 
field which is important only for long wavelengths and disappears in the high frequency 
limit, yielding again the null geodesic equation. 

3.1.3. Maxwell field: Let us now consider the Maxwell field with Lagrangian density 
and energy-momentum tensor 

Lm — ^TTT- Fc.d.F'^'^ , (27) 

"^"-b ~ ^ i^acFb'^ — ^ QabFcdF'^'^ ■ (28) 

The corrected conservation equation ([5]) yields 

(^FacFb' - i gabFcdF^"^ = - j^fj- F^cFb^V^R . (29) 

By taking now the high frequency limit 

Ab^Cbe''^ (30) 
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with the phase S{x'^) a rapidly varying function and C;, a slowly varying vector 
amplitude, and neglecting the derivatives of the latter, one obtains (using Fab — 

daA - dbAa), 



2aabS — Sa 







aab^'R 



where 5c = V^S", = CcC, = ScS", and 

Ctab — C^SaSb — {C^Sc) {SaCb + SbCa) 



S CaCb 



(31) 
(32) 

(33) 



2 



Sa — and 



Eqs. §^ and ^ yield 

and eqs. ([3T|) and ([32|) simplify to 
2aabS'' = , 

V''[(l + A/2)aa6] -0 

The only difference with respect to the standard Maxwell equations in curved space 
is the term in A in eq. ([36]). It is straightforward to see that, in this equation. 



(34) 

(35) 
(36) 



"afc 



A. 



» \aabV'' [l + A/2(i?)]| 



(37) 



where Xem is the wavelength of electromagnetic waves and L is the radius of curvature 
OS spacetime. In the high frequency limit << 1, the corrections to standard optics 
coming from the term A/2(-R) in the Lagrangian disappear and photons follow null 
geodesies and are transversal. In other words, the non-minimal coupling corrections 
to the Maxwell equations can only affect long wavelenghts, comparable to the radius 
of curvature of spacetime. In the high frequency limit, photons are transverse and 
propagate along null geodesies. 

3. 2. Perfect fluid with constant equation of state 

Let us now turn our attention to massive matter fields and, for simplicity, let us 
consider a perfect fluid composed of non-relativistic or relativistic particles with 
constant barotropic equation of state P = wp, where p and P are the energy density 
and pressure, respectively, described by the stress-energy tensor 

Tab (w + 1) PUaUb + Wpgab (38) 

and by the Lagrangian i„i = P = wp. The corrected conservation equation ([5]) yields 

{w + 1) UaU^VbP + (ui -I- 1) pU^VbUa 



+ [w + 1) UaV^Ub + WV aP = — 



\{W+\)F2 



pUaUb'S/ R 



l + A/2 

Projecting onto the direction of the fluid four-velocity u"" , one obtains 



^ + {w + l)pV''ub 



\{w + l) F2 DR 

1 + A/2 ^157 



(39) 



(40) 
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where r is the proper time along the timeUke fluid curves and D/Dt — u'^Vc- The 
correction to the conservation equation is best seen in the weak-field limit, in which 
T-t,^^§ = t+^/-Vp, and 



dp ^ r X{w + 1)F2 fdR ^ 



dt ' 1 + A/2 "^ydt 

Clearly, the fluid energy is not conserved and the right hand side of eq. ((1T|) acts as 
a source term describing the energy injected into the fluid per unit time and per unit 
volume. This correction disappears if A ^ 0, in vacuo, quantum vacuo {w — —1), or 
when DR/Dt = 0. 

We now project eq. ([39]) onto the direction normal to the four-velocity by the 
use of the projection operator h°'c defined by hat = gab + UaUh (recall that, in our 
signature, u'^Uc = — 1). This gives, after easy manipulations. 



Dt dr 
where 



r 



1 



Tt^u^u^ = r , (42) 
AFo 



-(i„-P)Vci? + VcP 



(43) 



.l + A/2 

The last term in square brackets, proportional to the pressure gradient, is the usual 
term that appears in GR and encapsulates the force exerted on a fluid element due 
to the fluid pressure (it is not to be attributed to the coupling between matter and 
R, nor does it signal any new effect). As mentioned earlier, a natural choice for the 
Lagrangian of a perfect fluid is Lm = P '30] . Substitution into eq. (|43|) yields 

„\ A ^cPW^, (44) 

i.e., remarkably, the only force is due to the pressure gradient already present in GR. 
For instance, in the case of dust with P = 0, it is = 0. Therefore, we conclude that 
there is no extra force which can be attributed to the coupling between the matter and 
R even for the case of a perfect fluid. The fact that energy is indeed not conserved for 
this fluid does not contradict this result. On the contrary, one can readily see from the 
right hand side of eq. ([55]) that the flow of energy only occurs along the direction of u"^, 
i.e., aligned with the fluid worldlines. Only the time component of the force is non- 
zero, while its spatial components in the frame comoving with the fluid always vanish. 
This kind of force can have no effect on the motion because, due to the normalization 
u'^Uc = — 1 , the only meaningful component of the four-acceleration a'^ and of the four- 
force cx a'^ is the one perpendicular to the four- velocity {a^u'^ = 0), i.e., the spatial 
one. For a timelike four-velocity u'^, a four-acceleration parallel to u'^ is necessarily 
zero, which is what we recover here. In a Newtonian analogy, the correction to the 
equations of motion would correspond to the introduction of a spatially homogeneous, 
but time-dependent, potential energy. 

The fact that no extra force appears for a perfect fluid was missed in Ref. |12j . 
We will discuss the implications of this fact, to some extent, in section [Jill For the 
moment, let us stress the following: In [T^] the authors use the opposite signature 
than the one used here, but also in Refs. [3U] where the result L,„ = P is derived. The 
use of a different signature is also the reason for the sign difference in the parenthesis 
and in front of P between eq. (|43p and the corresponding equation of Ref. 12]. Both 
equations are consistent given the signature adopted. However, when replacing L,„ in 
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this equation one has to take the signature into account as weU: since the gravitational 
action changes sign after a signature change, the sign of the matter action has to be 
changed for consistency, i.e., = —P with the signature convention of [H]. Note 
also that in a follow-up publication by Bertolami and Paramos studying the effects of 
the coupling discussed here on stellar configurations [33 , it is assumed that L„j = P, 
even though the same signature as in Ref. 12J is used, which is the opposite of the one 
used here and in Refs. [30]. This certainly affects the results derived and, therefore, 
one has to be cautious about the conclusions stated there until the effect of this 
inconsistency in the analysis is assessed. 

It is worth noting that the above conclusions apply to a massive scalar as well. 
This can be understood through the fact that a massive scalar admits a perfect fluid 
representation. Alternatively, it can be seen directly by the fact that the right hand 
side of eq. (O is proportional to Lmgab — Tab- Given the definition of the stress 
energy tensor, any part of the Lagrangian Lm which does not explicitly depend on 
the metric, such as the potential of a scalar field, does not contribute to this term. 
Massless or massive, a scalar field will always lead to an energy flow along its motion 
only, resulting into a vanishing extra force. 

3. 3. Conformal frames and energy conservation 

We have already discussed the fact that the i?-matter coupling can not be eliminated 
by a conformal transformation of the metric for a generic matter Lagrangian. This 
is evident in the fact that R couples in the same way to all the terms that a matter 
Lagrangian might be split into (kinetic, potential, etc.). However, one could think 
of employing the tool of conformal transformations in order to see how particular 
matter fields are affected by this coupling. More specifically, consider a matter field 
whose Lagrangian does indeed transform as L„i ^ ALm under the transformation 
gab — ^^dab, whcre A is some power of the scale factor il. Then it does seem possible, 
under certain conditions, to find a conformal frame in which energy is conserved for 
this field (and, consequently, particles associated with this field follow geodesies of 
this metric) . Recall also that geodesic motion is a characteristic of the theory, not of 
its representations, so changing conformal frames in an effort to find one that makes 
it manifest is perfectly legitimately Indeed, we can re-derive some of the results of 
the previous section with the use of conformal frames. This, besides being a way of 
verifying their validity, is also a nice exercise for realizing the ambiguities of the metric 
postulates, already pointed out in pi]. 

We begin by recalling that, if under the conformal transformation 

gab gab = ^^gab (45) 

the stress energy tensor Tab transforms as 

rpab ^ rpab ^srpab (46) 

where s is an appropriate conformal weight, then one can easily verify that the 
following equation relates the covariant divergencies of Tab and Tab .4 : 

Vaf"'' = n'VaT"'' + (s + 6)n'-^T''''Van - n''~^Tv''n , (47) 

* We are referring to whether a particle follows geodesies of some metric. The Einstein EP requires 
that the particles follow geodesies of some metric, not necessarily the one the action is written with. 
Of course, if particles corresponding to different fields follow geodesies of different metrics, then this 
signals a violation of the universality of free fall. Such a violation however, can only be detected 
in Eotvos-type experiments by comparing the trajectories of different particles and not by detecting 
deviation of a single species from geodesic motion (see section 13.41 1 . 
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where Vc is the covariant derivative associated with gab and T = gabT""^. Now, by 
imposing the condition that V aT"^^ — 0, i.e., that energy is conserved in the new 
frame, we get 

VaT"'' = ^ [Tg"'' - (s + 6)r''''] Vbfi . (48) 

We rewrite here eq. ([S]) for direct comparison 

A P 

^bT^' = Tv4t (^'"ff"' - ^"') • (49) 

i + A/2 

One can then easily notice that these two equations become identical if and only if 

s = -4, T = 2L,„, n^ = l + \f2{R). (50) 

The first two conditions give us characteristics of the matter fields and the last just 
pinpoints the appropriate choice of the conformal factor. The only case that we know 
of that satisfies the first two conditions (considering also the transformation rule (|46| 1 
is a perfect fluid with stiff equation of state P ~ p, a. particular realization of which 
is a massless scalar field (which admits a perfect fiuid representation). Therefore, 
if matter is only composed of such fluids or scalars, then a conformal frame can be 
found in which energy is conserved and the corresponding conformal factor is the one 
in eq. ((50)) . 

We stress once more that the fact that this procedure can not be carried out for 
all matter fields directly implies that energy conservation is not a generic characteristic 
of the theory. Also, the fact that massless scalar particles follow geodesies on the one 
hand, or the fact that we can not show geodesic motion for all fields on the other 
hand, should not mislead us to think that we can make any statement (positive or 
negative, respectively) about the EP. In fact, we have already shown above that also 
null dust, high-frequency photons, and perfect fiuid volume elements follow geodesies, 
even though the use of conformal frames could not provide these results. 

Since this is actually the lesson to be learned from this exercise, let us comment 
on it: In the perfect fluid (including null dust) case, we have already seen that energy 
is not conserved. Geodesic motion comes as a consequence of the fact that the energy 
flow is always timelike and parallel to the fluid worldlines instead of being purely 
spatial. For the case of the electromagnetic field, on the other hand, we already know 
that the matter action is conformally invariant, so one can not expect that much 
information can be obtained through the use of conformal frames. The bottom line 
is that reaching conclusions about geodesic motion and EP violations from energy 
conservation (the divergence of the stress energy tensor in some conformal frame) is 
not always an easy and unambiguous task as one may think when considering the 
metric postulates. This fact further supports the claims made in [21]. 

5.^. Equivalence principle violations and dark matter 

3.4-.1. Possible violations of the Equivalence Principle: As already mentioned, the 
unusual coupling between R and matter in the theories under investigation and the 
fact that energy is not conserved signal a violation of the EP. On the other hand, in 
order to get concrete evidence of that, but also in order to derive quantitative results, 
this issue should be examined further, since there are intricacies in using the metric 
postulates for judging whether, and how, a theory will show such violations. First of 
all, let us be more concrete and separate the Strong Equivalence Principle (SEP), the 
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Einstein Equivalence Principle (EEP), and the Weak Equivalence Principle (WEP) 
(see [3D] for definitions of these concepts). 

R couples in the same manner to all possible matter Lagrangians into which one 
could decide to split Lm (one for each matter field. Therefore, there is no reason to 
believe that test particles (which by definition do not contribute to the gravitational 
field) of different compositions will follow different trajectories. This implies that there 
is no strong indication that the WEP will be violated. However, there is an issue that 
should be approached with care: to which accuracy can a small particle be considered 
a test particle in these theories? Could it be that the coupling affects our ability to 
treat small particles as test particles and neglect their contribution in the gravitational 
field? This issue, which has received attention within the context of GR needs 
further investigation. 

On the other hand, there is not much to say at this point about the SEP. All 
theories with more degrees of freedom that GR violate the SEP, and the theory 
considered here is no exception, even without the i?-matter coupling. Obviously, 
constraints can be imposed by considering the evolution of the effective gravitational 
coupling. The i?-matter coupling should affect gravitational experiments, such as 
Cavendish experiments. An interesting point is that, since the coupling depends on 
the curvature, extended bodies whose size is comparable with the radius of curvature 
can get seriously affected. Considering the equation of motion of extended bodies in 
such theories will probably provide important constraints. 

Finally, one can consider the more accurately tested EEP. It would be interesting 
to see how the i?- matter coupling affects local non-gravitational physics. It should be 
stressed that this coupling can not be eliminated by a field redefinition, as mentioned 
above; therefore, going to a local frame can hold surprises. However, it is interesting 
to note that if f2{R) is such that /2(0) = then, assuming that spacetime is locally 
flat in the background by choosing suitable coordinates and treating matter as a 
perturbation, the effect of the coupling should make its presence felt only at second 
order and the relevant term should still be suppressed by A. 

All of the above are perspectives for future work. Let us close by discussing how 
the results derived here relate to that: We have found that massless particles follow 
null geodesies as usual. This essentially implies that experiments using light, such as 
redshift experiments, can not constrain the theory. We have also shown that perfect 
fluid volume elements follow geodesies of the metric g^^y. This includes dust, which is 
usually used as an approximation to study the motion of test particles. This result 
clashes with the claim that test particles will be affected by the i?-coupling, but the 
issue of how close a small particle is to a test particle remains open. Once more, it 
should be stressed that geodesic motion alone can not be used to draw conclusions 
about EP violations. However, one could say that our results hint towards the fact 
that violations of the EP might be more difficult to detect than expected. 

3.^.2. Is there an extra force that can account for dark matter? In [12j the theory 
under investigation was put forth as a resolution of the dark matter problem in 
galaxies. It was argued there that the extra force could account for the flat rotational 
curves of galaxies due to a "MONDian" behaviour. Even though in \l2i the extra 
force was calculated for a perfect fluid, its exact expression was not sufhciently 
detailed to support the argument of this MOND-like behaviour, which remained rather 
qualitative. However, our flndings cast doubts on whether the theory will really exhibit 
such a behaviour. As a matter of fact, if we model the galaxy as a perfect fluid (as 
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the authors of [H] do), then the extra force felt by the stars (the fluid-elements) 
vanishes and no dark matter phenomenology occurs. Another approximation (and 
a rough one as well) consists of considering a star in the outskirts of the galaxy as 
a test particle moving in vacuo under the influence of the gravitational field of the 
galaxy. In such an approximation, however, any correction to the motion with respect 
to the one predicted by Newtonian gravity would not be caused by some extra force, 
as claimed in jl2j : it would merely be due to the fact that the vacuum spacetime 
around the galaxy would differ from the one predicted by GR due to the difference 
in the field equations^ Therefore, it seems unlikely that an extra force due to the 
i?-matter coupling can actually explain the flat rotation curves of galaxies. 

4. Conclusions 

We have studied metric f{R) theories of gravity which include an additional direct 
coupling between the Ricci scalar and the matter Lagrangian. We have shown that 
they can not be cast into the form of either a usual scalar-tensor theory, i.e. one 
scalar field non-minimally coupled to gravity and minimally coupled to matter, or 
even a multi-scalar-tensor theory with unusual couplings between the scalar and the 
matter, such as those studied in |27l[26]. Therefore, they can not easily be dismissed 
as being equivalent to a theory with known phenomenology. On the contrary, we 
argued that the i?-matter coupling persists (at least implicitly as a coupling between 
R and a non-dynamical auxiliary field algebraically related to the matter) even after 
field redefinitions and can not be eliminated by the use of conformal transformations. 
Hence, the implication of its presence should be thoroughly examined. 

We took a first step towards such a study by examining possible deviations of 
the free fall trajectories from geodesies (of some conformal metric). Massless fields 
were considered, including a massless scalar field and the electromagnetic field, but 
in both cases motion followed null geodesies. The same result was obtained for null 
dust. Remarkably, also for a perfect fiuid or a massive scalar, the extra force vanishes 
and the motion remains geodesic as well. It is worth noting that this does not come 
as a consequence of energy conservation, but rather as a consequence of the fact that 
the energy flow is purely timelike and is aligned with the motion. 

These results seem to indicate that some of the simple tests of the EP might not be 
able to discriminate between the theory under consideration and theories that satisfy 
the EP. Phrased otherwise, this means that this theory might have better chances of 
escaping the tight constraints related with the EP than other EP-violating theories. 
However, the same results seem to cast doubts on whether this theory can actually 
achieve the goal for which it was initially put forth in Ref. [12], i.e., accounting for 
the observed galactic dynamics without resorting to dark matter. 

Although the last argument appears to discourage further study, our overall 
conclusions do motivate future work on the viability of theories with i?-matter 
coupling, since, besides the fact that they have already been used to address other 
problems [TH [71 [151 [S] as well, they also appear to be interesting toy models which 
can help us understanding issues related to the foundations of gravitation theory such 
as the EP. Additionally, the results presented here are, to some extent, preliminary: 

(t There have been attempts to address the dark matter problem in this way even in standard f{R) 
gravity without the R-matter couphng (see for example [35] )• Even though in vacuo the _R-matter 
coupling vanishes, the vacuum solutions of the theory considered here and those of f{R) gravity may 
differ due to the matching with an interior solution. 
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further work is needed to completely clarify the motion of massive and test particles in 
the theory, as well as its implications for Local Lorenz Invariance and Local Position 
Invariance. One of the main steps in this direction is the study of the post-Newtonian 
limit. Future work can also include stability issues, the Cauchy problem, as well as 
further investigations of the usefulness of the theory as a substitute for dark energy 
and dark matter in cosmology and astrophysics. 
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Appendix 

In this appendix we re-examine a recent claim by one of us in Ref. [28j that the 
action ([2]) could be cast into the form of a scalar-tensor theorjl^ by starting from the 
action 

S*^ y'd^xy^|M^ + i/((i?-</.) + [l + A/2(<^)]i™| , (51) 

instead of and attempting to show that it is also dynamically equivalent to 
When (|5ip is varied with respect to one obtains 

i(i?-0)/{' + A/^L,„ = O. (52) 

In the A = case (or if matter is absent), and provided that /{' ^ 0, this yields R = (j) 
which straightforwardly implies that for this specific choice for A (corresponding to 
the usual metric /(i?) gravity), the action (|5ip becomes dynamically equivalent to 
the reduced version of essentially the action ([1]). However, it was claimed in [55] 
that this equivalence continues to hold when A ^ (and matter is present). Clearly 
this can not be true, unless /2 — 0, in which case the i?-matter coupling ceases to 
exist. Indeed, in order for the action ([5T|) to be dynamically equivalent to eq. ([52]) 
should reduce to R = <j). This can not be achieved unless the second term is somehow 
proportional to R — (j) and the multiplying factor is strictly non- vanishing: this would 
require at best some unacceptable constraint on the matter Lagrangian, if possible at 
all. One can therefore conclude that the action (|5T|) can not be dynamically equivalent 
to ([2]), contrary to what is claimed in "28]. 

As a final note on field redefinitions with the use of an auxiliary scalar, let us just 
mention that, instead of starting from the action one could decide to begin from 

S" = yd4xy^|M^ + [l + A/2(i?)]L„ + *(i?-(/.)| . (53) 

Varying this action with respect to yields R — <j) and one recovers the action ^ as 
before. However, in this case, variation with respect to (j) yields 

* = In (54) 

II However, this claim was not actually used in order to obtain the main results of |28| . 



Modified gravity with R-matter couplings and (non-) geodesic motion 



15 



Assuming that /" ^ and replacing eq. ([7]) back in eq. yields 




0) 



(55) 



This, unlike the action ([5T|) . is dynamically equivalent to the action Using the 
field redefinition $ — /[{(/)) and introducing the potential V{^) = (</>/{ — /i)/2, the 
action ([55]) can take the form 



which is a Brans-Dicke theory with luq — (just like ordinary metric f{R) gravity) but 
with the addition of the unusual i?-niatter coupling. This way one can avoid having a 
non-linear gravitational Lagrangian at the price of introducing a scalar field without 
touching the i?-matter coupling in any way. 
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